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Abstract 

We find out the connection between the unintegrated gluon distribution at low intrinsic 
transverse momenta and the inclusive spectrum of the hadrons produced in pp collision at 
LHC energies in the mid-rapidity region and low hadron transverse momenta. The parame- 
ters of this distribution are found from the best description of the LHC data. Its application 
to the analysis of ep deep inelastic scattering allows us to obtain the results which describe 
reasonably well the HI and ZEUS data on the structure functions at low x. A connection 
between the soft processes at LHC and small x physics at HERA has been found. 



PACS number(s): 12.38.-t, 13.85.-t 



1 Introduction 



Hard processes involving incoming protons, such as deep-inelastic lepton-proton scat- 
tering (DIS), are described using the scale-dependent parton density functions. Usually, 
these quantities are calculated as a function of the Bjorken variable x and the square of the 
four-momentum transfer = —Q^ within the framework of the DGLAP evolution equa- 
tions [1] based on the standard coUinear QCD factorization. However, for semi-inclusive 
processes (such as inclusive jet production in DIS, electroweak boson production [2], etc.) 
at high energies which are sensitive to the details of the parton kinematics it is more ap- 
propriate to use the parton distributions unintegrated over the transverse momentum kt or, 
transverse momentum depend (TMD) distributions, in the framework of the /cj-factorization 
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QCD approach^ [3]. The /crfactorization formalism is based on the BFKL [5] or CCFM [6] 
evolution equations and provides a solid theoretical ground for the effects of initial gluon 
radiation. The unintcgrated gluon g{x,kt) (u.g.d.) and quark q{x,kt) distributions (u.q.d.) 
are widely discussed and applied in phenomenological calculations in the framework of the 
fct-factorization QCD approach and can be found, for example, in [7-23]^. In [2,10] the un- 
integrated parton distributions (u.p.d.) were obtained using the so-called KMR prescription 
within the leading order (LO) and next-to-leading order of QCD (NLO) at large Q"^ from 
the known (DGLAP-evolved [1]) parton densities determined from the global data analysis. 
These u.p.d. were successfully applied to analyze the DIS data at low x and a number of 
processes studied at the Tevatron and LHC (see, for example, [14-22]). However, at small 
values of Q'^ the nonperturbative effects should be included to evaluate these distributions. 
The nonperturbative effects can arise from the complex structure of the QCD vacuum. For 
example, within the instanton approach the very fast increase of the unintcgrated gluon 
distribution function at < kf < 0.5 GeV/c and = 1 (GeV/c)^ is obtained [11]. These 
results stimulated us to assume, that the u.g.d. in the proton can be determined also in the 
soft hadron production in pp collisions. 

In this paper we analyze inclusive spectra of the hadrons produced in pp collisions at 
LHC energies in the mid-rapidity region in a context including the possible creation of soft 
gluons in the proton. We estimate the u.g.d. function at low intrinsic transverse momenta 
kt < 1-5 — 1.6 GeV/c and extract its parameters from the best description of the LHC 
data at low transverse momenta pt of the produced hadrons. We also show that our u.g.d. 
similar to the u.g.d. obtained in [8] at large kt and different from it at low kt- The u.g.d. is 
directly related to the dipolc-miclcon cross section within the model proposed in [7, 8] (see 
also [9,28-32]) which is saturated at low Q or large transverse distances r ~ 1/Q between 
quark q and antiquark q in the qq dipole created from the splitting of the virtual photon 7* 
in the ep DIS. Here we find a new parametrization for this dipole-nucleon cross section, as 
a function of r, using the saturation behavior of the gluon density. 

The paper is organized as following. In Section 2 we study the inclusive spectra of hadrons 
in pp collisions and obtain the modified n.g.d. In Section 3 we discuss the connection between 
the u.g.d. and the dipole cross section. In Section 4 we apply the modified u.g.d. to describe 
the HERA data of the DIS structure functions: the longitudinal (Fl), charm {F2) and 
bottom (^2^) structure functions (SF). 

2 Inclusive spectra of hadrons in pp collisions 

2.1 Unintegrated gluon distributions 

The u.p.d. in a proton are a subject of intensive studies, and various approaches to inves- 
tigate these quantities have been proposed. At asymptotically large energies (or very small 
x) the theoretically correct description is given by the BFKL evolution equation [5] where the 
leading ln(l/a;) contributions arc taken into account in all orders. Another approach, valid 
for both small and large x, is given by the CCFM gluon evolution equation [6]. It introduces 

^See, for example, reviews [4] for more information. 

^The theoretical analysis of TMD was done recently in [24-27]. 
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angular ordering of emissions to treat the gluon coherence effects correctly. In the limit of 
asymptotic high energies, it is almost equivalent to BFKL [5], but also similar to the DGLAP 
evolution for large x ~ 1. The resulting u.g.d. depends on two scales, the additional scale 
g is a variable related to the maximum angle allowed in the emission and plays the role of 
the evolution scale /j, in the coUinear parton densities. In the two-scale u.p.d. obtained from 
the conventional ones using the Kimber-Martin-Ryskin (KMR) prescription [2, 10], the kt 
dependence in the unintegrated parton distributions enters only in the last step of evolution. 
Such a procedure is expected to include the main part of the collincar higher-order QCD 
corrections. Finally, a simple parametrization of the unintegrated gluon density was ob- 
tained within the color-dipole approach in [7, 8] on the assumption of saturation of the gluon 
density at low Q'^ which successfully described both inclusive and diffraction ep scattering. 
This gluon density xg{x, k^, Qq) is given by [8,9] 

xg{x, kt, Qo) = —^^^—Rl{x)k^exp (-Rl{x)k'^) , Ro{x) = (—) , (1) 

where ctq = 29.12 mb, = 0.2, Qo = 1 GeV, A = 0.277 and xq = 4.1 ■ 10"^ This simple 
expression corresponds to the Gaussian form for the effective dipolc cross section a{x, r) as a 
function of x and the relative transverse separation r of the qq pair [8]. In fact, this form can 
be more complicated. In this paper we study this point and try to find a parametrization for 
xg{x, kt, Qo), which is related to a{x, r), from the best description of the inclusive spectra of 
charge hadrons produced in pp collisions at LHC energies and mid-rapidity region. 

2.2 Quark-gluon string model (QGSM) including gluons 

The soft hadron production in pp collisions at not too large momentum transfer can be 
analyzed within the soft QCD models, namely, the quark-gluon string model (QGSM) [33-35] 
or the dual parton model (DPM) [36]. The cut n-pomeron graphs calculated within these 
models result in a reasonable description at small but nonzero rapidities. However, it has 
been shown recently [37, 38] that there are some difficulties in using the QGSM to analyze 
inclusive spectra in pp collisions in the mid-rapidity region and at the initial energies above 
the ISR one. However, it is due to the Abramovsky-Gribov-Kancheli cutting rules (AGK) [39] 
at mid-rapidity {y ~ 0), when only one-pomeron Mueller-Kancheli diagrams contribute to 
the inclusive spectrum pfi{y — 0,pt). To overcome these difficulties it was assumed [37] that 
there are soft gluons or the so called intrinsic gluons in the proton suggested in [40], which 
split into qq pairs and should vanish at the zero intrinsic transverse momentum {kt ~ 0) 
because at fcf ~ the conventional QGSM (without the intrinsic gluons) [33] is applied very 
well. The inclusive spectrum of hadrons Eda /d^'p (y ^ 0) = ph{y — 0,Pt) was split into two 
parts, the quark contribution pq{y ~ 0,pt) and the gluon one and their energy dependence 
was calculated in [37, 38] 

Ph{y - 0,Pt) = Pg{y ^ 0,pt) + pg{y ~ 0,pt), (2) 

where Pq{y ~ 0,pt) is the quark contribution and Pg{y ~ 0,Pt) is the contribution of gluons 
to the spectrum p^(y ~ 0,pt). It was shown [38] that the AGK cutting rules for the inclusive 
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spectrum at x ~ can be proofed within the QGSM [33] and pq{y ~ 0,Pt) can be presented 
in the form 

oo 

P.(O,Pt) = E^n(s)0^(O,Pt), (3) 

n=l 

where 0^ is the convolution of the quark distribution and the fragmentation function of the 
quark q to the hadron h that at x ~ (y ~ 0) is proportional to the pomeron number n, 
see (9) and (10) in [38]. Therefore, (3) is presented in the form 

oo 

Pg(0,Pt) = 4>q(Q,Pt) Yl ^(^n(s), (4) 
n=l 

where (T„ is the cross section for production of the n-pomeron chain (or 2n quark- ant iquark 
strings) decaying into hadrons, calculated within the "eikonal approximation" [41], the func- 
tion 0g(O,pt) is related to the pomeron-hadron vertex in the Mueller-Kancheli diagram. 
Inserting (7„ into (4) the quark contribution Pq{0,pt) of the spectrum is presented in the 
form [38] 

Pq{0,Pt) = 9{s/so)^(j)g{0,Pt), (5) 

where g = 21 mb, A = «p(0) — 1 ~ 0.12, ap(0) is the intercept of the sub-critical Pomeron, 
the function (/)q{0,Pt) is found from the fit of the SPS and LHC data on the inclusive spectra 
of charged hadrons in the mid-rapidity at the initial energies from 540 GeV till 7 TeV. 

The intrinsic gluons split into the sea qq pairs, therefore its contribution to the inclusive 
spectrum is due to the contribution of the sea quarks, which, according to the QGSM 
ideology [33], contribute to the n-pomeron shower at n > 2. Therefore, the contribution 
Pg(y ^ 0,pt) is presented in the form similar to (4) replacing 0^(0, pt) by 4>g{0,pt) and 
inputting n >2 

oo 

Pgi^^Pt) = 09(0, Pt) Y n(^n{s) = 
n=2 

/ ~ °° \ (6) 

M^^pt) [ Y ^^n(s) -Y(^n{s)\ = 

\n=l n=l / 

4>9{^,Pt){g{s/so)^ - and), 

where and = o"n(=^) is the nondiffractive cross section, the function (f)g(0,pt is also found 
from the fit of the SPS and LHC data on the inclusive spectra of charged hadrons in the 
mid-rapidity region [38]. The following parametrizations for 4>q{0,pt) and (f)g{0,pt) were 
found [37]: 

M^,Pt) = Aqexp{-bgpt), 
05(O,Pt) ^ Ag^tex.p{-bgPt), 

where sq — 1 GeV^, g' 21 mb, A = 0.12. The parameters are fixed [37] from the fit to the 
data on the pt distribution of charged particles at y = 0: Aq — 4.78 ± 0.16 (GcV/c)~^, bq — 
7.24±0.11 (GeV/c)-^ and Ag = 1.42±0.05 (GeV/c)-^; bg = 3.46±0.02 (GeV/c)-^ In Fig. 1 
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p, [GeV/c] 

Figure 1: The inclusive spectrum of the charged hadron as a function of pt (GeV/c) in the 
central rapidity region {y = 0) at a/s = 7 TeV at pt < 1.6 GeV/c compared with the CMS 
measurements [42] which are very close to the ATLAS data [43] . The dashed and dash-dotted 
curves correspond to the quark and gluon contribution given by (5) and (6), respectively. 
The solid curve represents their sum. 




p, [GeV/c] p, [GeV/c] 

Figure 2: The inclusive spectrum of charged hadron as a function of pt (GeV/c) in the 
central rapidity region {y = 0) at a/s = 540 GeV (left) and a/s = 900 GeV (right) compared 
with the UAl [44] and CMS [42] data. The long dashed curves are the quark contribution 
Pg{x = 0,pt) (5), the short dashed curves correspond to the gluon one Pg{x = 0,pt) (6), the 
dash-dotted curves are the sum of the quark and gluon contributions (2), the dotted curves 
correspond to the perturbative LO QCD [38]. The solid curves represent the sum of the 
calculations within the soft QCD including the gluon contribution (2) and the perturbative 
LO QCD. 

we illustrate the fit of the inclusive spectrum of charged hadrons produced in pp collisions 
at a/s = 7 TeV and the central rapidity region at the hadron transverse momenta pt < 1.6 
GeV/c. Here the solid line corresponds to the quark contribution pg, the dashed line is the 
gluon contribution pg, and the dotted curve is the sum of these contributions ph given by 
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Figure 3: The inclusive spectrum of charged hadron as a function of pt (GeV/ c) in the central 
rapidity region {y = 0) at = 2.36 TeV (left) and y/s = 7 TeV (right) compared with the 
CMS [42] and ATLAS [43] data. Notation of all curves is the same as in Fig. 2. 

(2). The little discrepancy between the data and our calculation (the dotted line) at pt > 1.2 
GeV/c disappears if the contribution of the perturbaive QCD within the LO is included. It 
is shown in Figs. 2 and 3, where the inclusive hadron spectrum is presented at ^/s = 0.54, 
0.9, 2.36 and 7 TeV, where the dash-dotted line is the sum of the spectrum ph{x = 0,pt) 
(2) and the short dash-dotted line is the result of calculation within the LO PQCD [38]. 
The quite satisfactory description of the data on such spectra was obtained in [38] at the 
SPS and LHC energies using (2) for ph{x = 0,pt) and the LO PQCD calculations both 
together. Therefore we conclude that the energy dependence of the inclusve spectrum of 
charged hadrons produced in the pp collision at the mid-rapidity region is reasonably well 
described using the parametrization (7) for 0^(0, pt) and (f)g{0,pt). 

2.3 Modified unintegrated gluon distributions 

As it can be seen in Figs. 1 — 3 the contribution to the inclusive spectrum at ?/ ~ due 
to the intrinsic gluons is sizable at low pt < 2 GeV/c, e.g., in the soft kinematical region. 
Therefore, we can estimate this contribution within the nonperturbative QCD model, similar 
to the QGSM [33]. We calculate the gluon contribution 0g(x ~ 0,pt) entering into (6) as 
the cut graph (Fig. 4, right) of the one-pomeron exchange in the gluon-gluon interaction 
(Fig. 4, left) using the splitting of the gluons into the qq pair. The right diagram of Fig. 4 
corresponds to the creation of two colorless strings between the quark/ ant iquark [q/q) and 
antiquark/quark {q/q)- Then, after their brake, qq are produced and fragmented into the 
hadron h. Actually, the calculation can be made in a way similar to the calculation of the sea 
quark contribution to the inclusive spectrum within the QGSM [33], e.g., the contribution 
(j)g{0,Pt) is presented as the sum of two convolution functions 

4>g{x,pt) = Fg{x+,Pht)Fq{X-,Pht) + Fq{x+,Pht)Fq{X-,Pht), (8) 

where the function Fq(^q){x+,pht) corresponds to the production of the hadron h from the 
decay of the upper vertex of qq string and Fg(^g){x-,pht) corresponds to the production of 



6 



h from the decay of the bottom vertex of qq string. They are calculated as the following 
convolution: 



Fq(g){x±,pht) = J ^ dxi J d^kit fg{g){xi, kit) Gq{q)^h (j^, Pht - kit)^ ■ 



(9) 




Cut 



=> 



c 

C 
C 

9 ^ 



Figure 4: The one-pomeron exchange graph between two gluons in the elastic pp scatter- 
ing (left) and the cut one-pomeron due to the creation of two colorless strings between 
quarks/ ant iquarks that decay into qq pairs, which are drawn as the semi-circles (right) [33]. 

Here Gg(^q)^h{z, kt) = zDg(^g)^h{z, kt), Dg(^g)^h{z, kt) is the fragmentation function (FF) of 



the quark (antiquark) to the hadron h, z = x±/xi, kf = Pht — h, x± = 0.5(ya;^ + x^ ± x), 

Xt = 2\J {m\ + Pt)/s. At X ~ 0, we get that x+ and x_ are equivalent to each other, e.g., 
x+ = x_ = mj/ a/s. The distribution of sea quarks (antiquark) /^(q) is related to the splitting 



function Vg^gq of gluons to qq by 

fq{q)i.Z, kt) = 



g{zi,kt,Qo)V, 



dzi 



(10) 

where 5'(-2i, ^it, Qo) is the u.g.d. The gluon splitting function Vg^gg was calculated within 
the Born approximation. In (10) we assumed the collinear splitting of the intrinsic gluon to 
the qq pair because values of kt are not zero but small. 

Calculating the diagram of Fig. 4 (right) by the use of (8) — (10) for the gluon contri- 
bution pg we took the FF to charged hadrons, pions, kaons, and pp pairs obtained within 
the QGSM [46]. From the best description of pg{x ~ 0,p/it), see its parametrization given 
by (7), we found the form for the xg{x, kt, Qo) which was fitted in the following way: 



xg{x,kt,Qo] 



3(Tn 



-Ci(l 



'X 



X 



(Rl{x)kl + C2{Mx)ktY) exp [-Ro{x)kt - d{Ro{x)kt 



where Rq{x) is defined in (1). The coefficient Ci was found from the following normalization: 

g{x,Ql)= 1"^° dklg{x,klQl), (12) 
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k, [GeV/c] 



Figure 5: The unintegrated gluon distribution xg{x, kt,Qo)/co (where cq = 3ao/ {4:TT^as{Qo))) 
as a function of kt aX x = xq and Qo = 1 GeV/c. The sohd and dashed curves correspond to 
the the modified u.g.d. (11) and original GBW gluon density [7] given by (1), respectively. 

and the parameters a = 0.7, C2 — 2.3, A = 0.22, hg = 12, d = 0.2, C3 = 0.3295 were found 
from the best fit of the LHC data on the inclusive spectrum of charged hadrons produced in 
pp collisions and in the mid-rapidity region as it can be seen in Figs. 1 — 3. 

In Fig. 5 we present the modified u.g.d. obtained by calculating the cut one-pomeron 
graph of Fig. 4 and the original GBW u.g.d. [7] as a function of the transverse gluon mo- 
mentum kt. One can see that the modified u.g.d. (the solid line in Fig. 5) is different from 
the original GBW gluon density [7] at small kt < 1.5 GeV/c and coincides with it at larger 
kt- This is due to the sizeble contribution of pg in (6) and (7) to the inclusive spectrum p{pt) 
of charged hadrons produced in pp collisions at LHC energies and in the mid-rapidity region 
(see the dashed line in Fig. 1). 

Let us also note that, as it was shown recently in [45], the modified GBW given by 
(11) describes reasonably well the HERA data on the proton longitudinal structure function 
Fi(Q2) at low X. 

3 Saturation dynamics 

According to [7,8] (see also [29-31]), the u.g.d. can be related to the cross section 
a{x, r) of the qq dipole with the nucleon. This dipole is created from the split of the virtual 
exchanged photon 7* to qq pair in ep deep inelastic scattering (DIS).The relation at the fixed 
Ql is the following [8]: 



Using the simple form for xg{x,kt) given by (1) as input to (13) on can get the following 




(13) 
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form for the dipole cross section: 



^GBw{x,r) = (To < 1 - exp 



(14) 



However, the modified u.g.d. given by (11) inputted to (13) resuhs in a more comphcated 
form for a(x, r): 



(Toil- exp 



bir 



Mix) 



(15) 



where bi = 0.045 and &2 = 0.3. In Fig. 6 we show the difference between the dipole cross 
section acBwi^ = xo,r) [8] and (Tinodif(a; = xo,r) obtained from the modified u.g.d. given 
by (11). The saturation effect means that the dipole cross section becomes constant when 
r > 2i?o- As it was shown in [7,8], at low the transverse 7*p cross section calculated 
within the dipole model is about constant when QRq < 1. It means that ~ I/-R0 can be 
treated as the saturation scale. One can see in Fig. 6 that the modified dipole cross section 
is saturated a little earlier than the GBW one. 



e> 1.0 




r[fm] 



Figure 6: The dipole cross section a/ao at function of r. The solid and dashed 

curves correspond to our calculations and calculation of [8], respectively. 

There are different forms of the dipole cross sections suggested in [28-30,48]. The dipole 
cross section can be presented in the general form [7]: 



a(x,r) = cTo5'(f^), 

where f = r/{2RQ{x)). The function (7(f^) can be written in the form [28] 

g{r^) = f2 log (1 + ^) , 

or in the form [48] 

l-exp|-f2log(^^ + e)|. 



(16) 



(17) 
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0.0 0.2 0.4 0.6 0.8 1.0 

r[fm] 



Figure 7: The dipole cross section a/ao at function of r. The sohd, dashed, 

dash-dotted and dotted curves correspond to our calculations, calculations of [8], calculations 
of [28] and [48], respectively. 



where saturation occurs for larger r. In Fig. 7 we illustrate the dipole cross sections (t/ctq 
at X = Xo which are saturated at r > 0.6 fm, obtained in [28,47,48]. They are compared 
with the results of our calculations (solid line) given by (15). The solid curve in Fig. 7 
corresponds to the modified u.g.d. given by (11), which allowed us to describe the LHC data 
on inclusive spectra of hadrons produced in the mid-rapidity region of pp collision at low pt- 
Therefore, the form of the dipole-nucleon cross sections presented in Fig. 7 can be verified 
by the description of the last LHC data on hadron spectra in soft kinematical region. 

Comparing the solid curve ("Modified a") and dashed curve ("GBW a") in Fig. 7 one 
can see that a^odifixjv) given by (15) is saturated earler than d'GBw{x,r) given by (14) 
with increasing the transverse dimension r of the qq dipole. If Rq = (1/GeV) (x/xq)'^''^, 
according to [7,8], then the saturation scale has the form Qg ~ I/-R0 = Qsoixo/^)^^"^ ^ where 
Qso = 1 GeV= 0.2 fm~^. The saturation of the dipole cross section (14) sets in when 
r ~ 2Rq or Qg ~ {Qsq/2){xo/x)'^^'^). Comparing the saturation properties of the modified a 
and GBW a presented in Fig. 7 one can get slightly larger value for Qso in comparison with 
g,o = 1 GeV. 



4 Proton structure functions 



The basic information on the internal structure of the proton can be extracted from the 
process of deep inelastic ep scattering. Its differential cross-section has the form: 



dxdy xQ^ 



l-y+y-\F2{x,Q' 



(19) 



where F2{x, Q^) and -Fl(x, Q^) are the transverse and longitudinal proton structure functions, 
X = Q'^/2{p ■ q) and y = Q^/xs are the usual Bjorken variables with p, q and s being 
the proton and photon four-momenta and total ep center-of-mass energy, respectively. In 
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the present paper we will concentrate on the charm and beauty part of F2{x,Q^) and on 
the longitudinal SF Fl{x,Q'^). Theoretical analysis [52-57] have generally confirmed that 
F2{x, Q"^) and F2{x, Q^) data can be described through perturbative generation of charm 
and beauty within QCD. The longitudinal SF Fl{x,Q'^) is very sensitive to QCD processes 
since it is directly connected to the gluon content of the proton. It is equal to zero in the 
parton model with spin 1/2 partons and has nonzero values in the framework of pQCD. 

In the /ct-factorization approach [3], the study of charm and beauty contributions to the 
proton SF F2{x, Q'^) and longitudinal SF Fl(x, Q"^) has been performed previously in [14-16]. 
In these calculations the different approaches to evaluate the unintcgratcd gluon density in 
a proton have been tested and a reasonable well agreement with the HERA data has been 
found, in particular, with the CCFM-evolved gluon density proposed in [9]. Below we apply 
the unintegrated gluon distribution given by (7) to describe recent experimental data [49,50, 
58-62] on F|(x,Q2)^ F^{x,Q'^) and taken by the HI and ZEUS collaborations at 

HERA. The main formulas have been obtained previously in [14]. Here we only recall some 
of them. 

According to the /cf-factorization prescription, the considered proton SF can be calculated 
as a following double convolution: 

F',\x,Q') ^Jjl de,C2{x/y,klQ',f,')f,{y,klt,'), (20) 

F^{x,Q') ^J2^}J^J dk^CL{x/y,klQ',f,')f,{y,klf,'), (21) 

where e'j is the electric charge of the quark of flavor /. The hard coefficient functions 
C2,l{x, k^, Q^, /i^) correspond to the quark-box diagram for the photon-gluon fusion subpro- 
cess and have been calculated in [14]. Numerically, here we set charm and beauty quark 
masses to rUc — 1.4 GeV and rub — 4.75 GeV and use the LO formula for the strong couphng 
constant as{i^'^) with Uf — A quark flavours at Aqcd = 200 MeV, such that as{Mz) — 0.1232. 
Note that in order to take into account the NLO corrections (which are important at low 
Q"^) in our numerical calculations wc apply the method proposed in [51]. Following [16,51], 
we use the shifted value of the renormalization scale nj^ = K Q^, where K = 127. As is was 
shown in [51], this shifted scale in the DGLAP approach at LO approximation leads to the 
results which are very close to the NLO ones. In the case of /ct-factorization this procedure 
gives us a possibility to take into account additional higher-twist and non-logorithmic NLO 
corrections [16]. 

The results of our calculations are presented in Figs. 8 — 10 in comparison with recent HI 
and ZEUS data [49,50,58-62]. Note that the data [49,50] on the longitudinal SF Fl{x,Q^) 
refer to the fixed value of the hadronic mass W — 276 GeV. The solid and dashed curves 
correspond to the results obtained using the modified u.g.d. (7) and the original GBW 
gluon, respectively. One can see that the predictions obtained with the u.g.d. (7) are in a 
reasonable agreement with the available HERA data for F2^{x, Q"^) as well as for longitudinal 
SF Fl{x, Q'^). Moreover, the shape of measured SF at moderate and large x values at high 

is better reproduced by the modified u.g.d. (7) compared to the original GBW one. When 
X becomes small, the predictions of both gluon distributions under consideration practically 
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Figure 8: The longitudinal structure function Fl{x,Q'^) at fixed W = 276 GeV and /x^ = 
K ■ Q^, where K = 127 [51]. Notation of all curves is the same as in Fig. 5. The HI data 
are taken from [49,50]. 
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Figure 9: The charm structure function F2c{x,Q'^). Notation of all curves is the same as in 
Fig. 5. The data are taken from [58-60,62]. 

coincide. Therefore we conclude that the link between soft processes at the LHC and low-a; 
physics at HERA is found, when we use the modified u.g.d. obtained from the description 
of pp-spectra at the LHC for the analyses of the behaviour of proton structure functions 
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Figure 10: The bottom structure function F2b{x,Q'^). Notation of all curves is the same as 
in Fig. 5. The data are taken from [58,59,61,62]. 

at HERA. Of course, it will be important for further studies of small-x physics at hadron 
colliders. 

5 Conclusion 

We have fitted experimental data on the inclusive spectra of charged particles produced 
in the central pp collisions at high energies taking into account the sum of the quark pg and 
the gluon contributions Pg in (2) and (7). The parameters of this fit do not depend on the 
initial energy in wide energy interval. Assuming creation of soft gluons in the proton at low 
transverse momenta kt and calculating the cut one-pomeron graph between two gluons in 
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colliding protons we found a form for the unintegrated gluon distribution (modified u.g.d) 
as a function of x and kt at fixed value of Qq. The parameters of this u.g.d. were found from 
the best description of the LHC data on inclusive spectra of charged hadrons produced in 
the mid-rapidity pp collisions at low pt. It was shown that the modified u.g.d. is different 
from the original GBW u.g.d. obtained in [8] at kt < 1.6 GeV/c and it coincides with the 
GBW u.g.d. at A;* > 1.6 GeV/c. 

Using the modified u.g.d. we have calculated the qq dipole-nucleon cross section CTmodif 
as a function of the transverse distance r between q and q in the dipole and have found 
that it saturates faster than (Tqew obtained within the GBW dipole model [7, 8] . Moreover, 
we have shown that the relation of the modified u.g.d. and amodif supports the form of the 
dipole-nucleon cross section and the property of saturation of the gluon density. 

It has been shown that the modified u.g.d. results in a reasonable description of the 
longitudinal structure function Fl{Q'^) at the fixed hadronic mass W . The calculations of the 
charm -^2^(2;, Q^) and the bottom -Ff (a;, Q"^) structure functions, within the /c^-factorisation, 
with using the modified u.g.d. and the GBW u.g.d. show a not large difference between the 
results for ^2^(2;, Q"^) in the whole region of x and and some difference for ^2^(2;, Q^) at 
low X and large Q^. The use of the modified u.g.d. results in a better description of the 
HERA data at large and small x. Therefore the link between soft processes at the LHC 
and low-x physics at HERA has been found, since the modification of the u.g.d. leads to a 
satisfactory description of both the LHC and HERA data. 
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